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Qs^ ' Abstract. In this paper we calculate the density profiles of virialized halos both in 

i^ , the case of structure evolving hierarchically from a scale-free Gaussian 5— field having a 

O ■ power spectrum P{k) oc fc" in a fi = 1 Universe and in the case of the CDM model, by 

Q using a modified version of Hoffman & Shaham's (1985) (hereafter HS) and Hoffman's 

H 
"^ (1988) model. We suppose that the initial density contrast profile around local maxima 



> 



is given by the mean peak profile introduced by Bardeen et al. (1986) (hereafter BBKS), 
and is not just proportional to the two-point correlation function, as assumed by HS. We 



X 

^ show that the density profiles, both for scale-free Universes and the CDM model, are not 



a 



power-laws but have a logarithmic slope that increases from the inner halo to its outer 
parts. Both scale- free, for n > —1, and CDM density profiles are well approximated by 
Navarro et al. (1995, 1996, 1997) profile. The radius a, at which the slope a = —2, is a 
function of the mass of the halo and in the scale-free models also of the spectral index n. 

Key words: Cosmology: theory - large-scale structure of Universe - Galaxies: formation 

1. Introduction 

The collapse of perturbations onto local density maxima of the primordial density field is 
likely to have played a key role in the formation of galaxies and clusters of galaxies. The 
problem of the collapse has been investigated from two points of view, namely that of the 
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statistical distribution of the formed objects (a question related to the biasing problem) 
(Kaiser 1984; Davis et al. 1985; BBKS) and that of the structure of these objects and 
its dependence on the statistical properties of the primordial density field (Gunn & Gott 
1972; Gunn 1977; Filmore & Goldreich 1984; Bcrtschinger 1985; West et al. 1987; HS; 
Hoffman 1988; Efstathiou et al. 1988; Quinn et al. 1986; Warren et al. 1991; White & 
Zaritsky 1992; Evrard et al. 1993; Crone et al. 1994; Navarro et al. 1995, 1996, 1997; 
Avila- Reese et al. 1998). 

To overcome the problem of the excessively steep density profiles, p ex r~^, obtained 
in numerical experiments of simple gravitational collapse Gunn & Gott (1972), Gott 
(1975) and Gunn (1977) were able to produce shallower profiles, p oc r^^ through the 
secondary infall process. Self-similar solutions were found by Fillmore & Goldreich (1984) 
and Bcrtschinger (1985), who found a profile of p oc r^^-^^. HS considered a scale-free 
initial perturbation spectra, P{k) oc fc" and assumed that local density cxtrema are the 
progenitors of cosmic structures and that the density contrast profile around maxima is 
proportional to the two-point correlation function. They thus showed that p oc r^" with 
a ~ {/Xnl ' ^^^^ recovering Bertschinger's (1985) profile for n = and = 1. They also 
showed that, in an open Universe, the slopes of the density profiles steepen with increasing 
values of n and with decreasing il, reaching a profile p oc r^^ for 17 ^ 0. Hoffman 
(1988) refined the calculations of HS and made a detailed comparison of the analytical 
predictions of the secondary infall model (hereafter SIM) with the simulations by Quinn 
et al. (1986) and Quinn & Zurek (1988). In spite of the high level of simplification of 
the SIM, and although the formation of dark matter halos in the numerical simulations 
seems to grow in mass by mergers with typically less massive halos, while the collapse 
in the SIM is spherical symmetric and gentle, those numerical simulations (Quinn et al. 
1986; Frenk et al. 1988) were in agreement with the predictions of the SIM. The good 
results given by the SIM in describing the formation of dark matter halos seem to be 
due to the fact that in energy space the collapse is ordered and gentle, differently from 
what seen in N-body simulations (Zaroubi et al. 1996). If things go really in this way, it 
is possible that galaxies and clusters of galaxies retain memory of their initial conditions. 
A great effort has been dedicated to study the role of initial conditions in shaping the 
final structure of the dark matter halos; but, if on large scales (evolution in the weakly 
non-linear regime) the growing mode of the initial density fluctuations can be recovered 
if the present velocity or density field is given (Peebles 1989; Nusser & Dekel 1992), on 
small scales shell crossing and virialization contribute to make the situation less clear. 
To study the problem, three-dimensional large-scale structure simulations were run with 
often confiicting results. While Quinn et al. (1986) and Efstathiou et al. (1988) found a 
connection between the density profiles of collapsed objects and the initial fluctuation 
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spectrum for Einstein-de Sitter universes [in particular Efstathiou et al. (1988) found 
density profiles steepening with increasing spectral index n], West et al. (1987) arrived 
at the opposite conclusion. In any case, the previous studies showed that the mass density 
profiles steepened with decreasing il. This result is in agreement with that of HS. More 
recent studies (Voglis et al. 1995; Zaroubi et al. 1996) showed a correlation between 
the profiles and the final structures. Finally Dubinski & Carlberg (1991), Lemson (1995), 
Cole & Lacey (1996), Navarro et al. (1996, 1997) and Moore et al. (1997) found that dark 
matter halos do not follow a power law but develop universal profile, a quite general profile 
for any scenario in which structures form due to hierarchical clustering, characterized by 
a slope /3 = -g^i^ = — 1 near the halo center and /3 = —3 at large radii. In that approach, 
density profiles can be fitted with a one parameter functional form: 
PJr) ^ Sn ^^ 

Pb Lh+Lf 

where ph is the background density and 5n is the overdensity [below we shall refer to 
Eq. (H) (Navarro et al. 1997) as the NEW profile]. The scale radius a, which defines 
the scale where the profile shape changes from slope /3<— 2to/3>— 2, and the 
characteristic overdensity, 5n, are related because the mean overdensity enclosed within 
the virial radius r^ is ~ 180. The scale radius and the central overdensity are directly 
related to the formation time of a given halo (Navarro et al. 1997). The power spectrum 
and the cosmological parameters only enter to determine the typical formation epoch 
of a halo of a given mass, and thereby the dependence of the characteristic radius on 
the total mass of the halo. Also these last results are not universally accepted. Recently, 
Klypin et al. (1997) and Nusser & Shet (1998) challenged the claim that a one parameter 
functional form could fit the density profiles using N-body simulations. Klypin et al. 
(1997), by using N-body simulations of CDM-like models, showed that the scatter about 
a one parameter fit is substantial, and that more than just one physical parameter is 
needed to describe the structure of halo density profiles in agreement with Nusser & 
Sheth's (1998) conclusions. In short, the question of whether galaxies and clusters mass 
density profiles retain information on the initial conditions and the evolutionary history 
that led to their formation remains an open question. 

In this paper, we introduce a modified version of HS and Hoffman's (1988) models to 
study the shapes of the density profiles that result from the gravitational collapse. In 
particular, we relax the hypothesis that the initial density profile is proportional to the 
two-point correlation function, and use the density profiles given by BBKS. 
The plan of the paper is the following: in Sect. 2 we show the reasons why the HS (1985) 
model must be improved and how it can be improved. In Sect. 3 we introduce our model 
and in Sect. 4 we show our results and finally in Sect. 5 we draw our conclusions. 
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2. Limits of the SIM and Hoffman & Shaham's approaches 

Most analytic work has focusscd on studying the evolution of isolated spherical systems 
(Gunn & Gott 1972, Gott 1975; Gunn 1977; FiUmore & Goldreich 1984; Bertschinger 
1985; HS; Hoffman 1988; White & Zaritsky 1992) because there is no analytical technique 
to study the evolution of a system starting from general initial conditions. It is known 
that, if the initial density profile of a halo is a power law in radius, S{r) oc r^™, {S{r) is 
the mean density at a distance r from a peak) then the density profile of the collapsed 
halo is also a power law: 

p oc r^" (2) 

where a = yS^, if m > 2, and a = —2 if m < 2 (Filmore & Goldreich 1984; Bertschinger 

1985). The restriction on m can be relaxed if non-radial orbits are permitted (White & 

Zaritsky 1992). HS and Hoffman (1988) related this spherical solution to dark matter 

halos, that form from initial gaussian density fields. In their paper they noted that the 

density profiles of very high peaks for r > ry have the same radial dependence as the 

correlation function of the initial field, (5 oc ^ ex r^(^+"). Then HS set m = 3 + n, 

obtaining: 

3(3 + n) .„. 

(4 + n) 

for n > —1, and a — —2 for n < —1. This last conclusion {a — —2 for n < —1) is not 

a direct consequence of the HS model, but it was an assumption made by the quoted 

authors following the study of self-similar gravitational collapse by Fillmore & Goldreich 

(1984). In fact, as reported by the same authors, in deriving the relation between the 

density at maximum expansion and the final one (see next section) HS assumed that each 

new shell that collapses can be considered as a a small perturbation to the gravitational 

field of the collapsed halo. This assumption breaks down for n < —1. In reality, all the 

slopes of the power spectrum satisfying the condition n < — 1 constitute a big problem 

for SIM (Zaroubi et al. 1996). In fact, as shown in Zaroubi's et al. (1996) simulations and 

Lokas' et al. (1996) paper, the slope n = — 1 marks the transition between two different 

dynamical regimes. In the cosmological context, such regimes correspond to a primordial 

perturbation field whose power spectrum is dominated by the high wavenumber modes, 

n > —1, or by the low wavenumber modes, n < — 1, respectively. 

In the regime n < — 1, the energy of the particles changes violently in time, and one 

expects at the end a state of statistical equilibrium, similar to the violent relaxation 

described by Lynden-Bell (1967). In this regime there is a strong dependence on the 

boundary conditions and the dynamics strongly depends on the last collapsing shell 

(Zaroubi et al. 1996). According to Lokas et al. (1996) for n < —1 the fluctuations grow 

faster than in the linear theory. 
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In the n > — 1 regime, the dynamics of the collapsed shells is hardly affected by the 
ongoing collapse of more distant shells. In this case, order is preserved in energy space 
and according to Lokas et al. (1996) a slowdown in the growth rate of perturbations is 
expected. In such a case, the SIM is expected to be a useful tool for calculating the final 
virialized structure of collapsed halos. For these reasons in this paper we suppose that 
-1 <n < 0. 

Let us also add that Eq. (||) is true only for regions outside the virial radius of a dark 
matter halo (see Peebles 1974; Peebles & Groth 1976; Davis & Peebles 1977; Bonometto & 
Lucchin 1978; Peebles 1980; Fry 1984). In the inner regions of the halo, scaling arguments 
plus the stability assumption tell us that ^{r) oc r c^+i) and we expect a slope different 
from Eq. (|3|). Syer & White (1996) found for the inner regions of the halo a profile 
p(r) (X r cs+i) , coincident with the slope of the correlation function. Nusser & Sheth 
(1998) found ^^±^ < a < ^^^, while Sheth & Jain (1996) found a = ji^- In 
other words, HS's (1985) solution applies only to the outer regions of collapsed halos and 
consequently the conclusion, obtained from that model, that dark matter halos density 
profiles can be approximated by power-laws on their overall radius range is not correct. 
It is then necessary to introduce a model that can make predictions also on the inner 
parts of halos. 

Another problem of HS's work is the assumption that S{r) ex ^(r) oc r~('^+"^. This is 
true only for very high peaks (see Ryden & Gunn 1987). As the peak height decreases, 
the peak profile becomes steeper than the correlation function (BBKS) and consequently 
the final density profile becomes steeper. Moreover, according to BBKS, the mean peak 
profile depends on a sum involving the initial correlation function, ^(r) oc 7-^(5+")^ and 
its Laplacian, \7^^{r) oc r"(3+«) (BBKS; Ryden & Gunn 1987). This means that there 
are at least two reasons why the density profile outside the virial radius must be steeper 
than in HS's model: 

— a) peaks that give origin to structures have height v = 2,3, and not ly —> oo; 

— b) S{r) depends on a sum of ^(r) and its Laplacian, and for peaks j/ = 2, 3 is steeper 
than the correlation function. 

Inside the virial radius HS's model cannot be used to predict the virialized density profile 
because it was a priori constructed to estimate density profiles at r > r^. 
That model suffers also from another drawback: namely, the assumption that the ac- 
creting matter is not clumpy, while in hierarchical scenarios for structure formation halo 
grows in mass by merger with less massive halos. As previously told according to Zaroubi 
et al. (1996) probably this last problem is not so difficult it, compared with the other 
two. 
In the following, we want to show that the predictive power of the SIM is greatly im- 



6 Del Popolo et al.: Density profiles of dark matter halos 

proved when the problems reported at point a) and b) are removed. These two problems 
are not problems of the SIM but only problems introduced by the HS's implementation 
of it. 

3. The model 

In the most promising cosmological scenarios, structure formation in the universe is 
generated through the growth and collapse of primeval density perturbations originated 
from quantum fluctuations (Guth & Pi 1982; Hawking 1982; Starobinsky 1982; BBKS) 
in an inflationary phase of early Universe. The growth in time of small perturbations 
is due to gravitational instability. The statistics of density fluctuations originated in 
the inflationary era are Gaussian, and can be expressed entirely in terms of the power 
spectrum of the density fluctuations: 

Pik) - {Kh (4) 

where 

<5k = / d^ kexp{-ikx.)S{x.) (5) 

^(x) = ^W^^ (6) 

Pb 

and ph is the mean background density. In biased structure formation theory it is assumed 

that cosmic structures of linear scale Rf form around the peaks of the density field, 6(x), 

smoothed on the same scale. 

If we suppose we are sitting on a i>a cxtremum in the the smoothed density field, we 

have that: 

,5(0) ^ J^e(O)'/' = va (7) 

together with: 

VS{r) |r=o- (8) 

If the Laplacian of S{r) is unspecified, that means that the extremum may be a maximum 
or a minimum, the mean density at a distance r from the peak is then: 

S{r) - J^e(r)/?(0)i/2 (9) 

(Peebles 1984; HS). If we calculate the mean density around maxima, as done by BBKS, 
by adding the constraint: 

\7^6{r) |,=o< (10) 

we find that the mean density around a peak is given by: 



{Sir)) 



vS,{r) '&iv-/,-f) 



^(0)1/2 ^(1-^2) 



2 



few + ^v^eM 



C(0)-i/2 (11) 
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(BBKS; Ryden & Gunn 1987), where v is the height of a density peak, ^(r) is the two- 
point correlation function: 
1 



i{r) 



P{k)ksiiii{kr)dk 



7 and _R» are two spectral parameters given respectively by: 
Jk^P{k)dk 



[Jk^P{k)dkJk(^P{k)dk] 



1/2 



R^ = 



3jk^P{k)dk 



1/2 



^k^P{k)dk 

while 19(71^, 7) is: 

3(1 - 72) + (1.216 - 0.974) exp [- (^) {^y 



(12) 

(13) 
(14) 



s{vin) = 



1/2 



(15) 






3(l-72) + 0.45+(^)^ 

In order to calculate 5{r) we need a power spectrum, P{k). In the following, we restrict 
our study to an Einstcin-De Sitter (f7 = 1) Universe with zero cosmological constant and 
scale-free density perturbation spectrum P{k) 



P{k) = Afe" 



(16) 



with a spectral index in the range — 1 < n < 0, and also to a CDM Universe with 
spectrum given by BBKS: 



P{k) ^ Ak-^ [In (1 + A.lQAk)f 



192.9 + 1340/c + 1.599 x Wk^ + 1.78 x 10'' fc^ + 3.995 x iC/c 



^6,4N-l/2 



(17) 



We normalized the spectrum by imposing that the mass variance at 8h~^Mpc is cts = 
0.63. In the case of a scale-free power spectrum, it is easy to show that the two-point 
correlation function can be expressed in terms of the confluent hypergeomctric function, 
iFi, and of the F function as: 

2X 1 



1 / ^ 



r(^-)iFi( — ,2'i^) 



where /3 — -R? /2, being Rf the filtering radius, and the Laplacian of ^(r) as: 



1 



2^2 -^ V 4/3; 2/3("+5)/2 



and finally 5{r) is 
<5(r) = 



n + 5 . -n-2 3 r2 

r(^-)iFi(^^,-,- 



(18) 



(19) 



6'(z^7,7)7 



1 



■ cxp 



1 



e(0)i/2 (1-72)^(0)1/2; 27r2 ^V 4/3; 2/3(»+3)/2 
n 3 r2\ e{iyj,j)Rl 



■Tl^\.,F, 



2 '2' 4/3; 37(1-72)^(0)1/2 



1 

2^ 



exp 



1 



4(3 J 2/3("+5)/2 



n + 5 



^1 



-n-2 3 



2' 4/3 



(20) 
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In the case that v is very large Eq. (15) reduces to 



^ 3(1 - 7')/(i.7) (21) 

and the mean density is well approximated by Eq. (0), which is the approximation used 
by HS to calculate (5(r). In reality, for peaks having i^ = 2, 3, 4, the mean expected density 
profile is different from a profile proportional to the correlation function both for galaxies 
and clusters of galaxies (see BBKS). For example for galaxies the CDM profile is steeper 
than that proportional to ^(r) as shown by Ryden & Gunn (1987) with a discrepancy 
increasing with decreasing v. As shown by Gunn & Gott (1972), a bound mass shell 
having initial comoving radius x will expand to a maximum radius: 

r„ = x/^(r) (22) 

where the mean fractional density excess inside the shell, as measured at current epoch 
io, assuming linear growth, can be calculated as: 

S^^rS{y)y'dy (23) 

'^ Jo 

At initial time ti and for a Universe with density parameter fli, a more general form of 

Eq. (H) (Peebles 1980) is : 

Tm = n = r (24) 

The last equation must be regarded as the main essence of the SIM. It tells us that the 
final time averaged radius of a given Lagrangian shell does scale with its initial radius. 
Expressing the scaling of the final radius, r, with the initial one by relating r to the turn 
around radius, r„j, it is possible to write: 

r = Frm (25) 

where .F is a costant that depends on a: 

F = F{a) = 0.186 + 0.156a + O.OlSa^ + QMla^ - 0.0045a'* + 0.0032a^ (26) 

(Zaroubi et al. 1996). If energy is conserved, then the shape of the density profile at 
maximum of expansion is conserved after the virialization, and is given by (Peebles 1980; 
HS; White & Zaritsky 1992): 

pM=P.(-)'§^ (27) 

V r / dr 

The density profile is a function of three parameters: the spectral index n, the density 
parameter fi, and the height of the density peak, u. In the limit v >> 1, the overdensity 
5{r) is proportional to the two-point correlation function and the density profile is a 
function of n and U, only, and then the expected profile is that by HS. 
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4. Results and discussion 

By using the model introduced in the previous section we have studied the density profiles 
of halos in scale- free universes with — 1 < n < 0, and for a CDM model characterized by 
a BBKS spectrum. As previously quoted, the chosen range of n is dictated by the limits 
of the SIM and by the values of n interesting in the cosmological context. The results of 
our calculations arc shown in Fig. 1 — 6. 

In Fig. 1, 2, 3 we have calculated the density profiles of halos in a scale- free universe with 
n = — 1, n = 0, and in a CDM model. In Fig. 1 we show the density profile for n = — 1. 




Log(r/Mpc) 



Fig. 1. Density profile for a scale- free spectrum with n = — 1 (solid line) and the NFW 
fit (dashed fine). The mass of the halo is ~ 2 x IO^'^Mq. 



The dashed line is the NFW fit, while the solid line represents the profile obtained from 
our model. The NFW profile fits well the density profile, except in the inner part. In 
the case n — 0, see Fig. 2, the NFW (dashed line) gives a good fit to the the density 
profile (solid line) also in the inner part of the density profile. The situation for the halo 
obtained from the CDM spectrum, smoothed on clusters scales Rf = Bh^^Mpc, (see 
Fig. 3) is similar to that of the case n — —1. This slope is in fact similar to that of the 
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standard CDM power spectrum on cluster scales. 

In Fig. 4 we plot the slope of the density profile for several values of n and v. The 
fundamental aim of the picture is to show how for small values of ly {ly = 2, 3), out from 
the inner region, the slope is larger than HF result. We began by finding the HS solution, 
that was recovered as expected in the limit v >> 1 or equivalently S(r) ~ '^C(^)/C(0)^ ■ 
This solution is represented by the short-dashed line which coincides with the result by 
HS, namely a — a^"/ , indicating an increase in the slope a with increasing n. Moreover 
the value of a is independent on the radius chosen to compute the slope which means 
that halos are described by pure power-laws, as described in the HF. Because of the rarity 




-1 

Log(r/Mpc) 



Fig. 2. As in Fig. 1, but now n = 0. The mass of the halo is ~ 0.5 x IO^^Mq 



of extremely high peaks, most galaxies and clusters will form from peaks of height 2 or 
3 (T (BBKS; Ryden & Gunn 1987): so we repeated the calculation of a for these values. 
If we choose a value of z^ = 3, the logarithmic slope of the density profile, calculated at 
Ih^^Mpc (we calculated the slope at a fixed distance because the density profiles are not 
power-laws) is steeper for all values of n (solid line) than that obtained by HS, and it is 
well approximated by Shet & Jain (1996) (dotted line), a = (h+n) ' obtained using stable 
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clustering and neglecting halo-halo correlations. At the same time the dependence of a 
on n is weaker than that shown by HS. We compared our result to that of Shet & Jain 
(1996) because their simple analytical formalism describes reasonably well NFW profiles 
on scales r > O.lr^ (at least till to Ih^^ Mpc). The choice of calculating the slope, a, at 
Ih^^ Mpc is suggested from the consideration that at that radius NFW profiles are well 
fit by Shet & Jain (1996) model. 

Obviously changing the radius at which the slope is calculted this reflects on the relation 
Qf-n, because the density profiles described by our model (excluding the case ly » I) 
are not power laws. For values larger than Ih^^ Mpc, the value of the slope is obviously 
larger. In any case, we have to remember that the comparison with Shet & Jain (1996) 
model is displayed only to show that density profiles have larger slopes than those found 
by HF. 
For i^ = 2 the slope is even steeper than the previous case (long-dashed line) and it is 




-1 

Log(r/Mpc) 



Fig. 3. Density profiles for a CDM spectrum smoothed on a scale Rf — 5h ^Mpc (solid 
hue) and the NFW fit (dashed hue). 



well approximated by Crone's et al. (1994) result, which is also consistent with the results 
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by Navarro et al. (1997) (sec their Fig. 13). As in the previous case the dependence of a 
on n is weaker with respect to that shown in HS. More massive halos j/ = 3 have flatter 
density profiles than less massive ones in agreement with Tornien et al. (1997). 
In Fig. 5 and Fig. 6 we plot a/vy, the variation of the ratio of the scale parameter, a 
and the virial radius, r^, versus M/M^ for a scale- free spectrum and a CDM spectrum, 
respectively. We remember that according to Navarro et. al (1996, 1997), a is linked to 
a dimensionless "concentration" parameter, c, by the relation a ~ — and the parameter 
c is linked to the characteristic density, S^^, by the relation: 

""~ln(l + c)-c/(l + c) ^'^^> 

The behaviour of the characteristic density of a halo, increasing towards lower masses 

in all the cosmological models, supports the idea that the M^-S^ relation is a direct 

result of the higher rcdshift of collapse of less massive systems. For scale-free models this 

implies 5^ oc M^^''^"^^/^ (same scaling relating M* and the mean cosmic density at a 

fixed redshift z). 

In the scale- free case, masses are normalized by the characteristic mass M*, which is 

defined at a time t as the linear mass on the scale currently reaching the non-linear 

regime: 

Air 
MM = -jRlPbit) (29) 

where the scale i?, is such that the linear density contrast on this scale is 5{R^,) = 1.69. 
Once known that the mass variance, ctm, for a power spectrum P{k) oc fc" is given by 
Cm oc _R^(3+n) .^-^^ remembering our normalization (TM{8h^^Mpc) — 0.63, the value of 
M* for n = — 1 results to be M* = 6 x 10^^ Mq. In the CDM case, the normalizazion 
mass, M*, is obtained imposing the condition croi-^-fl^)] = 1.69(1 + z). For our adopted 
normalization of the CDM power spectrum, M, = 3 x IO^^Mq, at z = 0. 
In the scale-free case. Fig. 5, the value of the scale radius a correlates strongly with halo 
mass and with spectral index n. The solid line represents a for n = —1. As shown in 
the figure, more massive halos have a larger scale radius a, or equivalently less massive 
halos are more concentrated. The dotted line shows a for n — 0. Also for this value of n 
more massive halos are less centrally concentrated. Finally from Fig. 5 we also see that in 
models with more small-scale power (or equivalently larger values of n) the haloes tend 
to have denser cores. These results were expected because halos with mass M << M, 
form much earlier than haloes with M >> M, and then arc more centrally concentrated. 
Moreover, for a fixed value of A'//M, , haloes form earlier in models with larger values 
of n and then have denser cores. This result is in qualitative agreement with those by 
Navarro et al. (1997), Cole & Lacey (1996), Tormcn et al. (1997). The filled squares and 
the filled exagons, in Fig. 5, represents a/r^ for n = —1 and n = respectively obtained 
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-0.8 



-0.6 



-0.4 



-0.2 



Fig. 4. The slope a of density profiles as a function of the spectral index n and v. The 
short-dashed line represents a in the limit z^ >> 1. ft coincides with the HS result. The 
solid line represents the logarithmic slope for v — i, while the dotted line is Shet & Jain's 
(1996) result. The long-dashed line represents a for v — 2. 

by Navarro et al. (1997) in the case / = 0.01 (see their paper for a definition of this 
parameter) which give the best fit to the results of their simulations. 
The virial radius r^ is obtained by using Navarro et al. (1997) equation: 



1.63 X 10" 



M 



1/3 



r!o 



-1/3 



{X + zY^h^^kpc 



(30) 



where flo is the actual value of the density parameter and h = 0.5. The virial radius 
determines the mass of the halo through: 

47r 



M, = 200pby r3 



(31) 



In the case n = — 1, our model gives less concentrated halos till M ~ lOM* and after 
this value the tendence is reversed. Maximum deviations of r-^ 2.5, between Navarro 
et al (1997) data and our model, are found in the low mass domain, < 0.25M*. The 
situation is similar to that described by Tormen et al. (1997) in his most re^aajed dynamical 
configuration (Fig. 16 second row in the left panel). In the case n = 0, our model gives 
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halos sligthly more concentrated in the overall studied mass range. In this case the 
discrepancies reach values of '-~^ 3. 
Also in the CDM case, Fig. 6, the value of the scale radius a correlates strongly with halo 




Log (M/M, 



Fig. 5. Trend of the scale radius a versus the mass of the halos in the case n = — 1 (solid 
line) and n = (dotted line). The filled squares and the filled exagons represents a/r^^ 
for n = — 1 and n = respectively, obtained by Navarro et al. (1997) in the case / = 0.01. 



mass. The solid line represents a for our model and the filled exagons the N-body results 
of Navarro et al. (1996,1997). As for scale-free spectrum, large halos are significantly 
less concentrated than small ones. This trend, a decreasing with decreasing M, can be 
explained in terms of the formation times of halos (Navarro et al. (1996,1997)). The 
dependence of concentration on mass is weak: there is only a change of a factor of ^ 4 
while M varies by 4 orders of magnitude. In this case the discrepancies between data 
and model is ^ 2. 

The result obtained is a remarkable improvement of the SIM being it able to reproduce 
almost all the prediction of N-body simulations with discrepancies of the same magnitude 
of those shown in Tormen et al. (1997) (this for the case n — —1) and surely much smaller 
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than those found in Cole & Lacey (1996). Our model is based on spherical simmetry, 



-0.8 
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cd 

*=* -1 2 - 

o 
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-1.6 
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Log (M/MJ 

Fig. 6. Trend of the scale radius a (solid line) versus the mass of the halos in the case 
of a CDM smoothed on a scale Rf = 5h^^Mpc and with normalization erg = 0.63. The 
filled exagons represents a/vv, obtained by Navarro et al. (1996,1997). 

and as we previously stressed, halos accretion does not happen in spherical shells but 
by aggregation of subclumps of matter which have already collapsed. In other words it 
seems that the halos structures does not depend crucially on hierarchical merging, in 
agreement with Huss et al. (1998). The SIM seems to have more predictive power than 
that till now conferred to it. 



5. Conclusions 

In this paper we have developed an improved version of the HS model to study the 
structure of the dark matter halos. We assumed that the initial density profile is given 
by the average profile given by BBKS and, solving the spherical collapse model, we 
obtained the virialized density profile. Our results can be summarized as follows: 
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— a) Differently from HS's (1985) model the density profiles are not power-laws but 
have a logarithmic slope that increase from the inner halo to its outer parts. In the 
outer parts of the halo, the density profiles are steeper than that found by HS and are 
consistent with p oc r~^, while in the inner part of the halo we find p oc~ r^^ . The 
analytic model proposed by Navarro et al. (1995) is a good fit to the halo profiles. 

— b) The radius, a, at which the slope equals —2 is a function of the mass of the halo 
and of the spectral index n. Lower mass halos are more centrally concentrated than 
the higher ones. For a given mass Af , halos having larger values of n have denser 
cores. 

— c) The good agreement of our (spherical simmetric) model with several N-body sim- 
ulations lead us to think, in agreement with Huss' et al. (1998) paper, that the role 
of merging in the formation of halos is not as crucial as generally believed. 
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